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■ABSTRACT 


-j  Laser-induced  periodic  pattern  formation  has  been  observed  ■ 

'on  a  variety  of  substances.  In  particular,  low-power  lasers  i  ^ 
-jhave  been  used’ to  deposit  a  pattern  on  a,  metal  surface.  For  a  -  ‘ 

-'relatively  smooth  surface  grating,  this  pattern  can  be  explained  ^ 

_jin  terms  of  a  perturbative  solution  of  Maxwell's  equations.  t'  r  ! 
— iHowever,  as  the-surface  grating '-is -^enhanced  by  this  initial  de-  i  //•''/  ' 
Jposition,  the  perturbation  solution  breaks  down.  An  alternate  /I  i 
_;non-perturbative  solution  of  Maxwell's  equations  for  such  rough  j 

jsur faces  is  considered  here.  Moreover,  other  possible  mechanisms 
that  may  assist  pattern  formation  are  discussed,  such  as  fleld- 
-jenhanced  evaporation  and  surface  migration. 


Oi 

A  V#  t  i  / 
i-'t  j  Spsjtpl 


jINTRODUCTION 


(IV.) 


^  Over  the  past  several  years,  a  number  of  researchers  [1-17]  have  j. 

-Observed  periodic  patterns  developing  on  the  surfaces  of  solids  that  were  j- 
-exposed  to  laser  radiation.  These  patterns  have  been  seen  on  various  metals,— 
-jinsulators ,  and  both  doped  and  pure  semiconductors.  Most  of  these  experl-  p 
jnents  were  carried  out  with  NdtYAG  or  ND: glass  lasers  tuned  to  wavelengths  p- 
Jbetween  1.0  pm  and  10.0  pm.  To  obtain  these  patterns,  laser  power  densities  - 
_[of  between  10  MW/cm^  and  1  GW/cm^  were  employed.  Above  this  range,  surface  p 
-jnelting  would  occur,  while  below  this  level,  no  discernible  surface  pattern  {- 
_ was  observed.  } — 


J  This  laser-induced  periodic  surface  structure  has  several  characteristics. 
-jFirst,  the  surface  pattern  consists  of  a  number  of  parallel  grooves  running  U 
-perpendicular  to  the  E  field  of  the  incident  laser.  Second,  Oron  and  [- 

JSorensen  [8]  demonstrated  that  the  underlying  surface  lattice  has  no  effect  |- 

_ !on  the  shape  of  the  laser-induced  surface  structure.  On  the  other  hand,  ■ — 

-jlsenor  [6]  has  shown  that  the  existence  of  scratches  on  the  surface  can  have 
a  significant  effect  on  the  final  surface  structure.  In  fact,  if  the  surface 
Is  randomly  scratched,  no  periodic  surface  structure  will  form.  Both  Brueck 
and  Ehrlich  [12]  and  van  Drlel  and  his  colleagues  [13]  have  also  observed  - 
-the  development  of  a  weaker  secondary  pattern  parallel  to  the  incident  E 
field.  With  this  secondary  structure,  the  surface  takes  on  a  scalloped  t- 

appearance.  Finally,  these  surface  patterns  produce  large  drops  in  the 
reflectivity  of  several  metals  [3-4].  • 


w  Most  of  the  work  so  far  has  considered  laser-induced  periodic  structures 
on  surfaces  that  already  had  a  certain  composition.  Needless  to  say,  this 
necessitates  large  laser  power  densities  (greater  than  10  MW/cm^)  in  order  to 
provide  sufficient  energy  for  the  surface  atoms  to  rearrange.  However, - 
-Brueck  and  Ehrlich  [12]  showed  that  only  a  modest  amount  of  laser  power  (less 
-than  10  W/cm2)  is  needed  to  establish  the  periodic  surface  structure  when  the 
constituent  atoms  are  deposited  from  the  gas  phase.  In  particular,  they  r 
showed  that  Cd,  Zn  and  A1  in  organometalllc  compounds  could  be  photochemical ly 
deposited  on  Si  or  SiO^  substrates  in  ordered  arrays.  [ 


To  understand  these  observations  and  resolve  any  difficulties,  one  must 
consider  how  this  laser-induced  periodic  surface  structure  is  produced. 

There  have  been  three  attempts  to  explain  these  effects.  First,  since  many 
metals  have  a  plasma  frequency  around  the  frequency  of  the  incident  laser 
.‘radiation,  surface  plasmons  will  be  excited  [11-12].  These  plasmons  will 

_ !  couple  with  the  electromagnetic  field  to  form  a  surface  polariton.  Since 

.  this  polariton  produces  maxima  and  minima  in  the  ^  field  of  the  surface, 
ij^he  surface  atoms  will  rearrange  themselves  to  minimize  their  energy  in  this 
..  E  field.  Second,  Maracas  et^  al^  [7]  have  observed  that  the  periodic  surface 
4 structure  can  be  considered  a  standing  wave.  This  wave  would  have  a  velocity 
Jvery  close  to  that  of  a  longitudinal  acoustic  phonon  of  the  substrate.  Thus,- 
_!  the  surface  pattern  may  well  be  a  phonon  excited  by  the  laser  that  was  - 

I "frozen"  in  place  by  the  cooling  of  the  lattice.  Finally,  van  Driel  and  his 
_i  associates  [15]  have  proposed  an  extensive  theory  based  on  the  small  inhomo-  - 
_  geneities  that  exist ; in  the  surface  layer.  This  initial  surface  roughness 
_  will  interact  with  the  incident  laser  beam  to  produce  a  dipole  moment  in  the 
_ surface  layer  (a  "radiation  remnant").  The  field  generated  by  this  dipole  >- 
_  layer  can  interferejwlth  the  refracted  beam  in  the  substrate  below  the  sur-  | — 
face.  This  interference  will  lead  to  inhomogeneous  energy  absorption  and  |- 
Jthus  the  redistribution  of  the  surface  atoms.-'  r 

-j  ^  ^ 

_j  In  the  following- section,  the  theory  of  the  surface  enhanced  E  field 
^produced  by  relatively  smooth  surfaces  and  the  resultant  pattern  formation 
_iwill  be  viewed. -olhe  inapplicability  of  this  approach  to  rough  gratings  will  — 
_|be  demonstrated.  To  overcome  this  limitation,  a  non-perturbative  solution  of 
JMaxwell's  equations  will  be  discussed.  Finally,  our  results  will  be  pre-  }- 
_jsentcd  along  with  suggested  improvements.  U 


_ THEORY  OF  SMOOTH  GRATINGS  !- 

I 

h 

Incident  radiation  converted  to  surface  plasmon  excitations  via  surface  - 
.roughness  has  been  investigated  extensively  in  the  past  decade  [18].  The  !- 
.classical  method  of  Rayleigh  [19]  has  often  been  used  to  solve  the  problem  J- 
-of  weak  scattering  [20].  Brueck  and  Ehrlich  [12]  also  adapted  this  formalism- 

_ to  describe  the  phenomena  of  laser-induced  pattern  formation.  Here,  we  will  ^ 

-briefly  outline  the  theory  of  Brueck  and  Ehrlich  and  discuss  its  validation  j- 
_and  limitation.  I- 

H 

Let  a  light  wave  of  frequency  u  represented  by  an  electric  field  •- 

•t  -imt 

-H  '  <'V 

“!  f 

jbe  Incident  on  the  metal  surface  at  a  normal  angle,  on  the  xz-plane.  If  we 

^express  the  surface  roughness  in  terms  of  a  Fourier  expansion. 


C(x,y) 


Tu  e^ 
1  8 


.  g  “  (g.,.g„).  r,,  -  (x,y). 


and  assume  that  each  amplitude  ug  is  snail  compared  to  the  wavelength  of  the 
incident  light,  then  taking  the  dielectric  function  as  a  step  function, 

e((i),r)  -  c(a))0(z  -  O  +  0(E  -  z).  (3) 


we  can.write  down  the  solutions  'of-  the  KaxwelL’s  equations  separately  for  the 
-two  media.  In  particular,  for  g^  >  oo^/c^and  z  >  ^,  we  have 

1  -  +1  +1  e'^^g^,  (4)^ 

i  r  sp  ' 

2  2  2  2  ■ 
where  kg  ■  g  -  w  /c  .  By  imposing  the  boundary  conditions,  l.e.,  the  local  ’ 


tangential  components  of  the  electric  and  magnetic  fields  have  to  be  con¬ 
tinuous  across  the  actual  sui^face  z  =  C(x,y),  we  can  express  the  reflected 
and  surface  plasmon  fields,  E  and  ^  respectively,  in  terms  of  E,  by 
writing  r  sp  i 


=  (1  +  ikoO 


and  assuming  g^  <<  1  In  the  continuity  equations.  If  we  consider  only  the 
p-wave  scattering,  i.e.,  the  light  polarized  in  the  x-directlon,  the  ■- 

expression  of  Brueck  and  Ehrlich  for  E  can  then  be  easily  obtained.  If  we' 
assume  a  linear  relationship  between  t^i  film  growth  rate  and  the  local  i. 

_ intensity,  L 


dF  *  ‘  •  .  . -  •  •  (6)“ 

at  nor  tt..  -•  t  ■ 

“then  the  suggested  equation  of  the  grating  depth  Ug  is  obtained  by  assuming  p 
~a  spatial  profile  of  T(x)  =  T  +rUCos(gx).  and  neglecting  the  |E  term.  p 

The  foregoing  theory  was.  developed  under  the  assumption  that  the  sur- 
“face  roughness  is  small,  g5<<  l.  However,'  as  more  metal  is  deposited,  J^he  p 
“roughness  increases.  .In. the  perturbative  theory  of  Brueck  and  Ehrlich,  E 
is  proportional  to  the  roughness  ii.  ' 'Therefore,  the ‘first  consequence  of  f 

“^increased  roughness  is  that  [e  l?^can  no  longer  be  neglected.  This  point  j 

may  already  be  i^each^d  when  tM  grating  is  only  a  few ‘layers  deep.  Moreover,’ 
“the  damping  of  E  can  no  longer  be  measured  from  the  mean  surface.  An  p 

“Increase  in  rbug^Hess  can  lead  to 'sufficient  enhancement  of  ?  that  localized 
“evaporation  of  metal  atoms  could  play  an  important  role.  Tranifer  of  large  p 
“amounts  of  energy  to  the  metal-carbon  bond,  such  as  in  the  deposition  of  Cd  p 
“on  Si  due  to  the  dissociation  of  gaseous  CdCCH^),*  ®ay  also  induce  sufficient" 
translational  motion  along  the  surface  so  that'^aasorption  does  not  occur  at  I 
“the  dissociation  site.  Such  dissociation  under  several  possibly  physisorbed 
“organometallic  layers  could  lead  to  trapping  of  the  organic  radical  in  the 
“Vicinity  of  the  surface.  This  "cage"  effect  could  lead  to  the  reformation  ' 
“of  the  organic  metal  bond.  Finally,  the  effect  of  surface  roughness  on  the 
“dielectric  function  might  be  Important.  Such  increases  of  surface  area 
could  substantially  alter  the  magnitude  of  this  function,  which  in  turn 
“could  enhance  the  rate  of  pattern  formation. 

“GENERALIZATION  TO  ROUGH  GRATINGS 

—  The  foregoing  considerations  lead  us  to  conclude  that  an  explanation  of 
“deep  pattern  formation  requires  a  non-perturbatlve  treatment  of  plasmon 
(formation,  and  at  the  present  time  a  numerical  solution  of  Maxwell’s  equations 
"seems  inevitable.  Recently  [21]  there  has  appeared  a  solution  of  the  problem' 
"“of  a  square-well  metallic  grating  in  an  applied  electromagnetic  field  ■' 

‘|(FIG.  lb).  Such  a  solution  is  restricted  to  explaining  the  most  qualitative 
'  ^features  of  sinusoidal  grating  formation.  We  have  established  a  general  ' 
formulation  of  the  problem  which  is  capable  in  principle  of  handling  gratings'^' 
of  any  shape  or  depth.  ^ 

In  the  H  polarization.  Maxwell’s  equations  can  be  written  as  [22]  [ 


+  c(x)k;Hy-  ^ 


()lnE(x)],' 


"where  H  is  the  component  of  the  magnetic  field  H  in  the  direction  perpendl- 
cular  tX  the  lattice  vector  of  the  grating  (FIG.  1).  For  square-well  T 

^gratings,  H  is  separable,  _ 


f 


Hy(x.y)  =  Z(z)X(x)  .  (8) 

and  Eq.  (7)  can  be  written  as  two  coupled  first-order  differential  equations: 


t  [e(x)k^"-A2]X  =  0. 
3x 


(9b)'.. 


.When  the  applied  radiation  is  incident  normal  to  the  grating,  it  is  found  !_ 
ithat  A  is  the  solution  of  the  equation  t, 

!  '  ,  not  L'c-  1  t(.'  i ■  ’3.)  L 

j  1  -  cos  (Bd/2)  cos  (ad/2)  +  l/2(ea/B  +  g/ (.ae))  sin(Bd/2)  sln(ad/2)=  0,  (10)_ 

J  •  ‘  h 

_ Vhere  d  is  th. period  of  the  grating,  e  is  the  dielectric  constant  of  the  1_ 

igrating,  a  =  (k  2  _  yy2)l/2  ^nd  B  =  (ek  ^  -  A^)^/^.  i 

“1  o  .  o  .  ^  r 

L-_-  dr:-.-!  v: '  tr.i-  •  t.i.-.  ...  [_ 

The  problem  to  be  considered  (FIG.  Ic)  is  that  of  a  multilayered  grating, 
lln  each  layer  the  grating  is;periodlc  and  the  .general  solution  to  Eq.  (9)  for 
the  n-th  layer  may  be  written  as 

“i  ■  K/.  (X)  ft?  )  .  . .  (1 1 V 


■V  itt!  o;  .. 


where  the  coefficients  A  and  B  are  to  be  determined  by  the  boundary  conditions 
between  layers,  and  A^^  is  the  £-th  solution  to  Eq.  (10)  for  the  n-th  layer. 

6y  applying  the  boundary  condition  of  continuity  of  H  and  ^  between  layers,  _ 
jwe  can  establish  the  recursion  relations  for  A  and 


X  <=  E  ,  +  F  , 

n  n,n-l  n-1  n,n-l  n-1 

I  =  C  ,  +  D  ,B  ,  , 

n  n,n-l  n-1  n,n-l  n-1 


(12ay 


(12b) 


“Inhere  the  matrices  C,  D,  E,  F  are  determined  by  the  periodic  structure  of 
-pach  layer.  Furthermore,  we  can  establish  boundary  conditions  between  the 
Top  two  layers  and  the  two  lowest  layers.  In  the  topmost  (infinite)  layer, 
-^he  eigensolutlon  is  written  as  [21] 


-IknZ 


2,1/2 


^IkjYnX  +  (1-Y„  )  z] 


-In  the  lowest  (infinite  layer)  we  have 


2,1/2. 


^IkolYn’'  -  (^  -rn  ) 


h 

(iA)r 


.  ,  nA 

where  y  ”  T" 
'n  d 


J  It  is  now  possible  to  establish 

i - 4:  -  aS  +  b^--  -  -  - 

J  ^ 

i 

j  ■  b2  +  a^, 

Vhere  D  is  a  unit  vector.  In  addition. 


(15a) 


(15b) 


(16) 


where  N  labels  the  penultimate  (finite)  layer,  and  the  matrix  L  is  determined 
by  the  parameters  which  characterize  that  layer.  However,  by  applying  the 
recursion  relations  of  Eq.  (12),  it  is  possible  to  establish  L 


K  ■  +  ■'2* 

B„  .  KjS  +  KjS 


(17a) 


(17b) 


Rearranging  (16)  and  (17)  gives 


(J^  -;LK2)‘ 


(LKj  -  Jj)D. 


(18)- 

l. 


.  fand  by  employing  Eq.  (12)  once  more  we  arrive  finally  at  a  matrix  equation 


Ifor  the  vectors 


and  B. 


M 

“  nS2 


n,n-l 


F  ,]■ 
n,n-l I 


n,n-l  n,n-l 


.'RESULTS  AND  DISCUSSION 


We  have  made  a  preliminary  calculation  for  a  square-well  silver  grating.'- 
..'The  field  intensity  |E(z,x)p  is  averaged  over  x  at  each  of  two  values  of  I. 

_ 2  (2  =  0  and  z  *  -h) ,  and  the  ratio  R  =  |  E(z  s  -h)  p/  |E(z  “  0)1”*  is  plotted  — 

as  a  function  of  the  well-depth  (roughness  parameter)  h  (FIG.  2).  R  is  thus  _ 
.'the  ratio  of  intensity  in  the  "peak"  region  of  the  grating  to  that  in  the 
j"well"  region.  For  very  shallow  gratings  (h  <  0.5  nm),or  about  two  mono- 
.jlayers)  it  is  seen  that  lE(z  =  0)l^  >  |E(z  =  -h)  |2  ,  which  supports  the  pertur-  - 
_,bation  description  of  Brueck  and  Ehrlich.  This  would  explain  how  a  particular 

_ ^grating  establishes  itself  out  of  the  "noise"  of  microscopic  roughness.  It  — 

'also  explains  why  the  establishment  of  such  a  grating  for  low-power  lasers 
_jrequires  the  photodissociation  of  an  organometallic  compound  and  does  not 
-(occur  in  the  presence  of  metal  vapor  alone.  However,  at  a  certain  depth  L 
J(around  0.5  nm)  | E(z  =  -h) P  begins  to  increase  very  rapidly.  It  seems  here  L 
_|that  the  grating  structure  will  only  continue  to  reinforce  itself  by  a 

_ mechanism  other  than  photodeposition.  Such  a  mechanism  might  be  photor  — 

Jevaporation  of  the  grating  itself,  whereby  the  trough  regions  are  "excavated"  - 
by  the  field.  Finally,  we  note  that  as  h  ■<■  «>,  |E(z  =  -h)  0.  This  is 

Jconsistent  with  an  upper  unit  to  the  growth  of  the  grating.  Work  is  currently 
underway  to  determine  whether  these  qualitative  trends  are  consistent  also  i 
with  slnu-soidal  gratings.  [- 
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Chapt.  1.  r 


FIG.  1. 

(a)  Sinusoidal  grating. 

The  hatched  area 
represents  the  metal.  L 

(b)  Square-well  grating  ' 

showing  a  separation  L 
into  three  layers,  one[— 
of  which  is  periodic  [ 
in  the  x-direction,  j.. 
and  two  of  which  are 
uniform.  [ 

(c)  Generalization  of  the  I 

square-well  grating  in 
which  there  are  three 
periodic  layers.  i- 
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— ■  FIG.  2.  Log-log  plot  of  the  ratio  R  =  |E (z  = -h) I  ^/|E (z  =  0) | ^  against  [- 

-i  the  grating  depth  h,  for  a  silver  grating  with  d  =  1050  nm  and  . 

-■  X  =  700  nm.  The  horizontal  line  corresponds  to  R  =  1.  Above 

-1  the  line  the  intensity  is  greatest  in  the  region  of  the  well  ]- 

(z  =  -h) ;  below  the  line  the  intensity  is  greatest  in  the  p 

region  of  the  peaks  (z  =  0)  of  the  grating.  Clearly  the  latter- 
condition  pertains  for  very  smooth  gratings  (h  -»■  0)  and  also  !- 
for  very  rough  gratings  (h  -♦•  »)  [  - 
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